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ABSTRACT 

Emission from X-ray binaries (XRBs) is known to be a major component of the total X-ray luminos- 
ity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and 
evolution of X-ray binaries on very long timescales (~ 10 Gyr). In this paper, we present results from 
large-scale population synthesis models of binary populations in galaxies from z = to 20. We use 
as input into our modeling the Millennium II Cosmological Simulation and the updated semi-analytic 
galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history (SFH) 
and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the parameters 
known from previous studies to affect the evolution of XRBs. We use the results from each model, 
along with observationally derived prescriptions for hot gas emission, to calculate the integrated X- 
ray luminosity of each galaxy in the catalog and create galaxy X-ray luminosity functions (XLFs) 
for several redshift bins. We compare our models with observed galaxy XLFs from Tzanavaris & 
Georgantopoulos (2008) and find that some of our models are able to reproduce their overall shape, 
normalization, and evolution. We use the results from our highest likelihood model to track the evo- 
lution of the normal galaxy X-ray luminosity density out to z — 4 and find that its evolution is driven 
largely by XRBs in galaxies with X-ray luminosities between 10^° and W^^erg s~^. 
Subject headings: stars: binaries: close, stars: evolution. X-rays: binaries, galaxies, diffuse background, 
galaxies: stellar content 



1. INTRODUCTION 

X-ray binaries (XRBs) are believed to be major con- 
tributors to the overall X-ray luminosity of normal galax- 
ies (those not dominat ed by the emission of a nuclear su- 
perrnassive blac k hole ; iFabbiandllQSl iKim fc Fabbianol 
120031 : lFabbianoll2006 0. Normal early- type galaxies have 
older stellar populations and their X-ray emission is dom- 
inated by low-mass XRBs (LMXBs) and hot interstellar 
medium (ISM) . On the other hand, the X-ray emission of 
normal late- type galaxies, which are still actively forming 
stars, have significant contributions from both LMXBs 
and high mass XRBs (HMXBs). 

X-ray and multiwavelength studies of galaxies using 
Chandra and XMM-Newton have yielded a great deal 
of information about the X-ray luminosities of galaxies, 
including many X-ray correlations that ha ve been estab- 
lished to hold out to at least z = 1 (e.g. iLehmer et al.l 
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120071 120 081: iSvmeonidis et al.l I201H rVattakunnel et al.l 
[2012; Bas u-Zvch et all 120121: iCowie et al.ll2012D . These 
relations include a strong correlation between X-ray 
emission from HMXBs and the star forma ti on rate (SFR) 
of g a laxies (e.g. IGilfanov et al.l l2004t iLehmer et al.l 
IMot IMineo et al.l 120121 ) as well as a scaling relation 
between the emis sion from LMXBs and the stellar 
mass of a gal axv (IGilfanov eFall [2001 iGilfanovl [200l 
Leh mer et all 120101; iBoroson et al.l 120111: iZhang et al.l 
l2012f ). Recent ultradeep Chandra and multiwavelength 
surveys (e.g [Bra ndt & Hasingcr 2005) have allowed for 
robust tests of these relations in ver y distant galaxies. 
For example, iBasu-Zych et al.l (120121) use a 4 Ms expo- 
sure of CDF South (IXue e7ani201in and X -ray stacking 
to study faint X-ray sources out to z ~ 8, finding that 
the relation between X-ray production and star forma- 
tion rate undergoes a small amount of evolution out to 
z ~ 4 that is likely driven by the metalicity evolution of 
HMXBs. 

Galaxy XLFs derived from rece nt observat ions show 
significant evolution with redshift (iNorman et al. 200'3; 
Ptak et al.' '2007; 'Tzanavaris fc Geor gantopouloa .2008) . 
Tzanavaris & Georgantopoulos (2003) (T&G08 here- 
after) use data from three Chandra deep fields (CDFS, 
E-CDFS, and CDFN) and the wide area survey XBootes 
to compile observations of 207 X-ray luminous normal 
galaxies (101 early-type and 106 late-type) out to z ^ 
1.4. They find a clear evolution of the galaxy XLF nor- 
malization with redshift that is driven almost exclusively 
by late-type galaxies. More specifically, this evolution is 
proportional to (1 -I- z)'^ with fc = 2.2 ± 0.3 for the total 
population, fc = 2.4^21) for late-type galaxies, and, for 
the early-type population, fc = 0.7j]i[ g (consistent with 
zero). Because XRBs are major contributors to the total 
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X-ray emission of normal galaxies, observationally de- 
rived XLFs can put constraints on theoretical models of 
XRB formation and evolution. 

At present, there has been little theoretical work done 
on the evolutio n of XR B populations over cosmolog- 
ical t i mescales dWhite fc Gho sh 1 991 iGhosh fc Whit^ 
I2OOII: iZuo fc Lill20TT[ r it is thought that XRBs 
could play a major role in t he evolution of these 
XLFs (IWhite fc GhoshI [1991 IGhosh fc Whitd [200I 
INorman et all 120041 ). IGhosh fc White! (|200in . using a 
semi-empirical, semi-analytical approach, linked XRB 
lifetimes with star formation rates (SFRs), showing that 
SFRs that are evolving on cosmological timescales sig- 
nificantly affect the XRB populations and, therefore, the 
integrated galactic X-ray emission. This predicted evo- 
lution should be evid ent even at lower redshifts (z < 1; 
IWhite fc Gh"osbl[l998h . 

Recently, the advances in available multi-wavelength 
observations of distant galaxies, as well as our under- 
standing of binary stellar evolution and galaxy forma- 
tion and evolution, have reached a level of maturity that 
allows us to conduct an in-depth study of the XRB pop- 
ulations of distant galaxies. In this paper, we study the 
evolution of X-ray binaries on cosmologically significant 
timescales, using data from detailed, large sc ale simula- 
tions. We use data from a catalog created by iGuo et aTl 
(|2011f) using semi-analytical galaxy evolution models ap- 
plied to the recent Millennium Cosmological Simulation. 
These data are used in tandem with the binary popu- 
lation synthesis (PS) code, StarTrack, to simulate the 
XRB populations of individual galaxies from 2: = to 
~ 20, taking into account the full star formation and 
merger histories of each galaxy. From these models, we 
derive the integrated X-ray emission of each galaxy and 
compare the resulting galaxy XLFs and their evolution 
to observed galaxy XLFs. Our goal is to obtain better 
constraints on the parameter space for models of XRB 
formation and evolution, and to better understand the 
nature of t he X-ray emission of galaxies at high redshifts. 

Recentlv. iFragos et al. (2012) (F12 hereafter) used sim- 
ilar techniques as those described in this paper to study 
the evolution of the global XRB population with redshift. 
They model how the total Universal specific X-ray lumi- 
nosities (Lx per unit stellar mass, star formation rate, 
and volume) of LMXBs and HMXBs evolve over cosmic 
time out to z ^ 20. Their models were constrained by 
observed luminosities of HMXB and LMXB populations 
in the local Universe. They found that the LMXB pop- 
ulation dominates the total population at low redshifts, 
with HMXB contributions becoming dominant for red- 
shifts higher than z ~ 3.1. 

The outline of the paper is as follows. Section 2 de- 
scribes our simulation tools, StarTrack and The Millen- 
nium Simulation II, and the methodology we follow in 
developing our models of XRB populations in galaxies. 
Section 3 describes how we compare our models to ob- 
servational results, namely those of TfcGOS, and the sta- 
tistical analysis we use to determine our best models. In 
section 4, we describe and discuss our results, and we 
conclude with a summary in section 5. 

2. SIMULATING X-RAY LUMINOSITIES OF GALAXIES 
2.1. The Millennium Cosmological Simulation 



The Millennium Cosmological Simulation is an un- 
precedented computational effort to sim ulate the dark 
matter distribution in the Universe (see iSpringel et al.l 
120051 for details). In this study, we use the data from 
the most recent M illennium Run II (MRU hereafter; 
iBovlan-Kolchin et al. .2009) . This is an N-body simu- 
lation that follows the evolution of 2160^ particles each 
of mass 6.9 x 10^h~^Msun within a co-moving box with 
sides each of size 100h~^Mpc. The cosmological model 
used in the simulation is a ACDM model with il„i = 0.25, 
= 0.75, rib = 0.045, and h = iJo/100 km/s/Mpc = 
0.73. 

The MRU has 60 snapshots in time that were saved 
and analysis was done to identify substructures within 
the dar k matter di stribu tions, including dark matter 
halos. Guo et al.; (201l|, Gil hereafter) use a semi- 
analytic procedure to track the evolution of the galaxies 
that exist within these halos. Once subhalos are iden- 
tified, their merger trees are derived. The evolution of 
these subhalos provide the base for the galaxy formation 
m odel. The models u s ed by Gil build upon the work 
of iDe Lucia fc Blaizoti (|2007[) . making improvements to 
the treatment of supernova feedback, reincorporation of 
ejected gas, galaxy sizes, the distinction between satellite 
and central galaxies, and the effect of the environment 
on galaxies. 

The result of Gil's semi-analytic model is a catalog 
of the galaxy population at 60 different times between 
2 20 (about 13.4 billion years ago) and the present 
day. These catalogs include properties such as metallic- 
ity, stellar mass, bulge mass, the mass of hot and cold 
gas, rest frame luminosity magnitudes, etc., for each of 
the galaxies in the simulation box, as a function of time. 

2.2. StarTrack 

To simulate the XRB populations of the galaxies from 
MRU, we use StarTrack, a current binary population 
synthesis code that has been tested and calibrated us- 
ing detailed binary star calculations and incorporates all 
the most important physical proc esses of binary evolu- 
tion (|Belczvnski et al.l [2001 120081 ): 

(i) The evolution of single stars and non-interacting bi- 
nary components, from their birth, taken as the time 
of their initial emergence onto the main sequence, to 
compa ct remnant fo r matio n using evolutionary formu- 
lae of Hurlcv et al. (2000*) modified as described in 
iBelczv nski et al. (2008). Various wind mass loss rates 
and their effect on stellar evolution are also incorpo- 
rated into the code a nd have been recently updated 
(IBelczvnski et al.ll2010D 

(ii) The time evolution of orbital properties. StarTrack 
numerically integrates a set of four differential equations 
describing the evolution of orbital separation, eccentric- 
ity, and component spins, taking into account tidal in- 
teractions, magnetic braking, gravitational radiation and 
stellar wind mass losses. 

(iii) All types of mass-transfer phases. This includes both 
stable and unstable mass-transfer processes, which are 
driven by either nuclear evolution or angular momentum 
loss. Unstable mass transfer is encountered most often 
as a direct consequence of rapid stellar expansion during 
nuclear evolution, but angular momentum loss may also 
lead to instability. 

(iv) Supernova explosions, which are treated by taking 
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into account mass/angular momentum losses as well as 
supernova asymmetries (through natal kicks to neutron 
stars and black holes). 

(v) X-ray emission, which is tracked for accreting bina- 
ries with compact object primaries (both for wind-fed 
and Roche- lobe overflowing systems). The resulting X- 
ray luminosities are calculated from the secular averaged 
mass accretion rate, but are not calculated for unstable 
accretion phases because the timescales are very short. 

The models in this paper include a recent revi- 
sion of StarTrack that incorporates updated stellar 
winds and their re-calib rated dependence on metallicity 
(|Belczvnski et al.]|201ClD . However, two more recent up- 
grades have not been incorporated into these results, as 
the simulations were run long before the changes were 
made. The first update includes a revised neutron star 
and BH mass distribution based on fully co nsistent su- 
perno va simulations (BeJczynski et al, 2 0TTI : iFrver et al.l 
I2OI2I) . The second, most recent upgrade improves upon 
the treatment of donor stars in CE events via usage of 
the actual value of the A parameter, the measure of the 
central concentration of the donor and envelope binding 
energy, for which us ually a constant value is assumed 
(|Dominik et al.ll2012t ). 

Table 1 lists the input parameters of our population 
synthesis models, which can be put into two categories. 
In the first category there are the parameters that cor- 
respond to the initial properties of the binary popula- 
tion. These values are relatively well constrained by the 
most recent observational surveys. Also in the group are 
stellar wind prescriptions and natal kick distributions, 
which can also be constrained by observations. In the 
second group are the truly "free" parameters that cor- 
respond to poorly understood physical processes, which 
we are not able to model in detail. One of these truly 
"free" parameters is the common envelope (CE) effi- 
ciency (acs), which measures how efficiently orbital en- 
ergy loss is transformed into thermal energy that will 
expel the donor's envelope during the CE phase. We note 
that in our calculations we combine acE and A, the bind- 
ing energy parameter described above, into one CE pa- 
rameter. Whenever we mention the CE efficiency acE, 
we refer in practice to the product qq e x A, effectively 
treating acE x A as a free parameter (see lBelczvnski et al.l 
[20081 for details). 

We create a grid of 192 PS models, a subset of those 
used in F12, run for nine different metallicities and each 
simulating 5.12 x 10^ stars for 14 Gyr. In this grid, we 
varied all the parameters known from previous studies 
to affect the evolution of XRBs and the formation of 
compact ob jects in general ("Belczvnski et al." 2007ll2010l : 
iFragos et a l. 2008, 2010; Linden et al. 2009). Specifi- 
cally, we vary the CE efficiency, initial binary mass ratio 
distribution, initial mass function (IMF), supernova (SN) 
kicks for direct collapse BHs, and stellar wind strength. 
We also take into account the possibility of CE inspi- 
rals with Hertzsprung gap donors that could terminate 
binary evolutio n barring the subsequent XRB formation 
(jBelczvnski et al . 2007). 

For all models, we assume a Maxwellian distri bution 
of supernova kicks given by iHobbs et al] (|2005D . with 
<T — 265fcm/s. For compact objects formed with partial 
mass fallback, the natal kicks given by the IHobbs et al.l 
p005.) distribution are decreased by a factor of (1 — //b). 



where //& is the fraction of the stellar envelope that 
falls back after the SN explosion. In our standard pre- 
scription, direct collapse (DC) BHs, BHs formed with 
ffb — 1, are given no natal kick. However, due to re- 
cent theoretical evidence that even the most massive stel- 
lar bl ack holes have probably received small as ymmetric 
kicks iLinden et al.ll2010HVarsecchi et al.ll2010[ ). m some 
models used in this work we set a lower limit (0.1) on 
the amount by which the natal kicks may be decreased 
due to mass fallback, allowing for small natal kicks to be 
given to direct collapse BHs. 

The mass of the primary star in each binary is de- 
termined by the adopted IMF. It is important to note 
here that, because we sample the IMF with only the pri- 
mary star, we are only sampling the high mass end of the 
IMF because the primary stars that form XRBs must be 
massive enough to form a BH or NS. The mass of the 
secondary star is calculated using a distribution function 
for the binary mass ratio, q = M secondary I M primary- We 
vary the distribution of q between a flat distribution in 
the range q = to 1 and a distribution that has 50% of 
the binaries follow a distribution with g = to 1 and the 
other half follow a "twins" distribution, with q = 0.9 to 
1. 

The models and their numbers used in this work are 
the same as those used in F12, except here we exclude 
the models 97-192, which have all s ystems f ollowi ng the 
pure "twins" q-distribution. Frago s~t al.l ()2012[ ) show 
that models with the pure "twins" distribution are un- 
physical, as they fail to reproduce local populations of 
XRBs. Thus, we exclude them here. 

We want to note that our population synthesis code 
calculates the bolometric luminosity of each XRB based 
on the rate of mass transfer. In order to compare our 
model results with observed datasets we need to estimate 
the X-ray luminosity in a specific energy band, which 
in this study is the soft X-ray band of Chandra (0.5 — 
2.0 keV). In order to calculate the bolometric correction, 
we used two sets of published X-ray spectra from Galactic 
neutron star and black hole XRBs at different spectral 
states (iRemiUard fc McClintocki I200a IWu et al.l 12010) . 
Following the same procedure outlined in F12, we derive 
to the bolometric correction factors for different types of 
XRBs and use these results to estimate the 0.5 — 2.0 keV 
X-ray luminosity of our modeled XRB population. 

We also note that the PS models used here take into 
account only binary systems formed in the field and 
not those formed via dynamical interactions in dense 
stellar systems. Dynamical formation in globular clus- 
ters (GC), for example, is a significant formation path- 
way for LMXBs in old, GC rich elliptical ga laxies (e.g. 
iHumphrev fc Buotd [2Ciol iZhang et all 120121 ). Further, 
the LMXB populations formed in GCs can be as much 
as 2 — 3 times more luminous than the field population 
in bright elliptical galaxies (?). 

For each model and each metallicity value, we calculate 
the X-ray luminosity as a function of age for single bursts 
of star-formation, where we also take in to account the 
effect of transient XRBs (see iFragos et al . 2008, 2009). 
Taking into account the assumed initial mass function 
and initial mass ratio distribution, we normalize the total 
X-ray luminosity to a nominal population of 10^" Mq. 
This quantity, Lx,spec{t)i^rg/s/10^°MQ), is the specific 
X-ray luminosity of a single age stellar population as a 
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function of its age. The specific X-ray luminosity coming 
from our population synthesis models can be convolved 
with the star-formation history and metallicity evolution 
of a galaxy to calculate the total X-ray luminosity of its 
complex stellar population. 

2.3. Convolving StarTrack with the Gil catalog 

The MRU catalog created by Gil corresponds to 60 
snapshots that span a redshift range from z = to 
~ 20. For each galaxy, we can derive its complete pro- 
genitor tree. Each progenitor galaxy has a unique SFR 
and metallicity, so for every stellar population in a tar- 
get galaxy we know during what timeframe and at what 
metallicity that population was created. We then con- 
volve the star formation histories with Lx,spec{t) derived 
from the PS models for the appropriate metallicity val- 
ues. 

The SFRs given for each galaxy in the Gil catalog 
are averaged over the entire timestep. At, between sub- 
sequent snapshots so that the total new stellar mass cre- 
ated in a given progenitor galaxy is Mnew = SFRprog^t. 
In order to account for the possibility of starbursts, we 
assume that all new stellar mass forms in a 20 Myr burst 
occurring at a random time between subsequent snap- 
shots, ti < tburst < ^i+i, where ti is the timestamp as- 
sociated with the snapshot of a given progenitor galaxy. 
This effect is important only for the HMXBs of young 
populations, and the effect on LMXBs is minimal since 
their evolution occurs on timescales much longer than 
the timesteps between snapshots. A 20 Myr duration 
is reasonable, given that the most cited values for star- 
burst durations are around 3-10 Myr fe.g. lThornlev et"al] 
120001 : [Harris et al.. ,2004 ), though there is evidence for 
longer burst duration s on the order of a few 10* yrs (e.g. 
iMcQuinn et"aI|[2QTol ). 

By summing the soft band X-ray luminosities of all 
the stellar populations in a given galaxy, we derive the 
integrated X-ray luminosity from XRBs. The end result 
is a catalog of integrated XRB luminosities of galaxies 
within the MRU co- moving volume from z = to z ^ 20. 

3. COMPARING WITH OBSERVATIONS 

3.1. Galaxy Classification 

When comparing our results to the observations of 
T&G08, we want to distinguish between early and late- 
type galaxies. For their sample, T&G08 cross correlate 
with other surveys to obtain optical counterparts for 
their X-ray selected galaxies, which they used for clas- 
sification. 

The classification of a galaxy as early or late-type 
can be based either on its morphology or its spectro- 
scopic properties. In many o bservational surveys o f dis- 
tant galaxies (e.g. GOODS; iGiavalisco et al.l[200l ). the 
morphologies of most galaxies that are observed cannot 
be determined due to inadequate spatial resolution, and 
therefore colors are used instead. 

For color classification of our model galaxies as early or 
late-ty pe, we adopt the method developed by Bell et al.l 
(|2004l ) . They showed that it is possible to define the pop- 
ulation of early-type galaxies empirically by using the 
bimodality of the color distribution, which they stud- 
ied out to z 1. The MRU database includes absolute 
rest frame magnitudes in the SDSS ugr filters, which 



can easily be tr ansformed to the UBV filters used in 
IBell et alJ (l200l . The magnitudes include the effects of 
dust extinction. Following the IBell et al.l ()2004[ ) prescrip- 
tion, we define early type galaxies to be galaxies where 
{U -V) > 1.15 - 0.31z - 0.08 * (My - 5logioh + 20). 
Figure 1 shows plots of {U — V) vs. My — 5logioh for 
the MRU galaxies at various redshifts with the cutoff 
function overlaid on top and the bimodality is clearly 
present. 

If we instead define galaxy types based on morphology 
with late type galaxies having "^^"'"^ < 0.7 and early- 
type galaxies ^j^j^'J"^ > 0.7, we find that there is approx- 
imately a 1% and 0.5% contamination among the color 
defined late and early- types respectively. Thus, these 
methods give nearly identical results, but using colors to 
define morphology allows us to better simulate observa- 
tions. The morphology method of classification, since it 
is independent of color, provides a check on our color clas- 
sification. The fact that they both yield similar results 
is encouraging and indicates that the colors provided by 
the Gil catalog are able to yield accurate morphology 
classifications. 

3.2. Creating Model X-ray Luminosity Functions 

For our analysis, we only select galaxies with stellar 
masses greater than 10^ solar masses, as galaxies with 
mass less than this are very unlikely to have X-ray lumi- 
nosities that are observationally relevant. For instance, 
the dwarf galaxies in the SINGS sample, with masses 
~ IO^Mq, generally have X-ray luminosities below 10"^^ 
ergs / s (0.5 - 8 keV) with many having no binaries de- 
tected at aU above 5 x 10^^ ergs / s (0.5 - 8 keV) (Jenkins, 
et. al. 2012 [in prep]). The Gil catalog contains several 
hundred of these more massive galaxies at very high red- 
shift (z — 19.9) and over 2 million at z = 0. It should 
also be noted that we assume that all of the galaxies in 
our sample are normal galaxies. This is a valid assump- 
tion because galaxies with AGN only constitute ~ 2 — 4% 
of all galaxies and, there fore, any selection effects on our 
data would be minimal (iXue et a l. 2010; Haggard et al.l 
IMotlSilverman et al.ll2009D . 

To compare our results with observations, we derive 
the XLF for our galaxies by calculating the number den- 
sity of galaxies versus their integrated X-ray luminosity. 
The simulation data are all within a single co-moving 
volume of constant size. For the time slice represented 
in each snapshot, 0(L) is defined as: 

VmriiSl 

Here, Lmin,Lmax are the bin limits, Vmrh is the vol- 
ume of the MRU simulation, which is (100 Mpc/h)^, 
and Sl is the size of the luminosity bin in log space, i.e. 
Sl = LogiQ{jf^). Herein lies a major difference be- 
tween the theoretical and observational luminosity func- 
tions. Observational surveys study a range of redshifts 
within a light cone. An entire volume of space cannot 
be observed at a constant redshift so a range of red- 
shifts is explored. Thus, when calculating (/!)(i), ob- 
serve rs such as TfcG08 u se methods like the one found 
in iPage fc Carreral (|2000[ ) which uses the following defi- 
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Table 1 

Model Parameters 



Parameter 


Notation 


Value 


Reference 


Initial Orbital Period distribution 


F(P) 


flat in logP 


Abt (1983) 


Initial Eccentricity Distribution 


F(e) 


Thermal F(e) ~ e 


Heggie Q975) 


Binary Fraction 


/bin 


50% 




Magnetic Braking 






Ivanova & Taam f2003) 


Metallicity 


z 


0.0001, 0.0002, 0.005, 0.001, 






0.002, 0.005, 0.01, 0.02, 0.03 




IMF (slope) 




-2.35 or -2.7 


Kroupa (2001); Kroupa & Weidner (2003) 


Initial Mass Ratio Distribution 


F(q) 


Flat, twin, or 50% flat plus 50% twin 


Kobulnickv & Frver (2007): Pinsonncault & Stanek (2006) 


CE Efficiency 


"CE 


0.1, 0.2, 0.3, or 0.5 


Podsiadlowski ct al. (2003) 


Stellar wind strength 


^wind 


0.25, 1.0, or 2.0 


~ Bclczvnski et al. (2010) 


CE during HG 




Yes or No 


Belczvnski et al. (2007) 


SN kick for ECS/AICt NS 




20% of normal NS kicks 


Linden et al. (2009) 


SN kick for direct collapse BH 




Yes or No 


Fraeos et al. (^010^ 



^ Electron Capture Supernova / Accretion Induced Collapse 




-22 -20 -18 -16 -22 -20 -18 -16 -22 -20 -18 -16 
- 5 log, oh - 5 log,oh - 5 logioh 



Figure 1. The rest-frame U — V color of simulated galaxies against the absolute magnitude in V-band, My — blogioh. It is clear that the 
Gil galaxy catalog exhibits the same bimodality as observed galaxies. The red line corresponds to 1.15 — 0.312 — 0.08* {My — 5logiQh + 20). 
This may be compared to the red sequence fitting of U-V colors by Bell et al. (2004); see their figure 1. Galaxies that lie above the red line 
are considered to be early-type galaxies and those that lie below are late-type galaxies. 
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nition: 



N 



JL„i„ Jz. 



%dzdL 



(2) 



Here, ^ represents the rate of change of the survey vol- 
ume with respect to redshift and ZmaxiL), ZmmiL) are 
the redshift ranges for a source as a function of himinos- 
ity such that it stays within the flux limits of the survey 
and within the redshift interval. Our simulated galaxies, 
on the other hand, exist within a co-moving volume and 
our snapshots capture all galaxies that exist at a given 
redshift. In order to compensate for this difference, we 
adopt similar redshift intervals as used by T&G08 and 
calculate the XLF for each redshift individually using 
equation (1) above. This gives us 4>{L,z). Then, for 
each luminosity bin centered at L, we take the average 
value of </)(L, z) over all z in the interval, giving us an 
estimate for (^(L) for that bin in that redshift interval. 

3.3. X-ray Luminosity from Hot Gas 

In addition to XRBs, the hot ISM in a galaxy can have 
a significant contribution to its overall X-ray luminosity. 
T&G08 do not distinguish between emission from XRBs 
and hot gas in their analysis of X-ray bright galaxies. 
Their analysis is done in the soft X-ray band, so emission 
from the hot ISM becomes important and needs to be 
taken into account when calculating the total integrated 
luminosities of our galaxies. 

We use observationally derived relations to estimate 
the X-ray luminosity of hot gas in early-type galaxies 
from their K-band luminosity. It has been shown that 
there is a power law relationship between the X-ray lumi- 
nosity of the hot ISM in early-type galaxies and both the 
K-band luminosity (Lx oc L]:^^^) of the galaxy and the 
temp erature of its hot gas (T cx L^^^; iBoroson et al.l 
1201 if ). The K-band luminosities of MRU early- type 
galaxies are easily calculated from mass and age using 
synthetic stellar population models (Bcrtclli ct al. 200^. 
With these relations, we estimate both the full band X- 
ray luminosity and the temperature of the hot gas in 
each early-type galaxy. We then input the temperature 
into the APEC XSPEC X-ray spectral fitting package to 
calculate the k-corrected X-ray luminosity of the hot gas 
in the soft X-ray band. 

For late-type galaxies, we estimate the hot gas X-ray 
luminosity based on the SFRs given in the Gil cata- 
log and the power-law relationship catalog between the 
soft-band X-ray luminosity of the hot I SM and the SFR 
for lat e-type galaxies (Lx oc SFR^-^'^; [Strickland et al.l 
[200413). 

In summary, we estimate the total X-ray emission from 
hot gas in all of the galaxies in the Gil catalog. We add 
those values to each galaxy's integrated X-ray emission 
from XRBs, calculated using StarTrack, to obtain the to- 
tal integrated X-ray luminosity of each galaxy. We find 
that on average XRB emission contributes to 50 — 60% 
of the 0.5 — 2 keV emission from bright (Lx > lO^^^) late- 
type galaxies and ~ 40% of the 0.5 — 2 keV emission 
from bright early-type galaxies for our best fitting model 
(205). Hence, we find that hot gas emission has an ap- 
preciable effect on the galaxy XLFs. See section 4 for 
more details on our results. 



3.4. Statistical Analysis 

From deep Chandra survey observations, T&G08 
present early- and late-type galaxy counts in several lu- 
minosity bins and two redshift intervals, < z < 0.4 and 
0.4 < z < 1.4. Associated with each count is a survey vol- 
ume that depends on both luminosity and redshift. Let 



the set of counts he d = {di. 



1. 



.NJ = 1...M} 



and associated volumes be Vi,j, where i ranges over the 
N = 5 luminosity bins and j ranges over the M = 2 
redshift bins. T&G08 assume that the dij are subject to 
Poisson statistical errors. 

Similarly, for a particular choice of parameters, 0, (for 
example, see Table[2]) our model produces a set of counts 
n = {ui j \i = 1 . . . N,j = 1 . . . M} for galaxies in each 
luminosity and redshift bin in our (lOOMpc/ft.)^ model 
volume. We assume that the counts n are drawn from a 
Poisson distribution with (unknown) means A — {K.j\i — 
1 . . . N,j = 1 . . . M} (note: this A is separate from the 
one used before to describe the CE efficiency parameter) . 
Because we only observe one particular set of counts, n, 
we do not measure the rates A implied by our model di- 
rectly, but instead must treat A as a nuisance parameter 
whose distribution under the observed n must be inte- 
grated over. 

Bayes' rule relates the posterior probability of model 
parameters 9, p{9\d), to the likelihood of the data under 
the model, p(d\9), the prior probability of the model pa- 
rameters before the data have been observed, p{9), and 
a normalizing constant, p{d), called the evidence, that is 
independent of 9 via 



Pi9\d) 



Pid\9)pi9) 
p{d) ■ 



(3) 



Writing the likelihood in terms of the (unknown) true 
mean A implied by the model, we have 

P{d\d) - n/ dK,jp{d.,j\K,j)p{K,j\9) , (4) 



where 

Pid^,J\Xi,j) = ''"''^^''^^ exp(-WijAij) , (5) 

is the Poisson probability of drawing di,j counts in a vol- 
ume Vij = Vij (100Mpc//i)"^ when the underlying rate is 

Xij per (100Mpc//^)^ 

The distribution of the underlying rates implied by our 
model, p{Xij\9) must be estimated from the observed 
Uij. Applying Bayes' rule again, we have 

P(A^,^■|6') =p(A,j|n,j(6')) = X —■ (6) 

P ['T'iJ ) 

The counts observed in a model with underlying rate A^ ^ 
are Poisson distributed, so 



P in^,J\Xi,j) = — ^ exp (-A,j) . 



(7) 
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Table 2 

Parameters and likelihood values for models referred to in this paper. A full list is available online. 



Model 




IMF exponent 


'Vwind 


CE-HG 


q distribution 


DC BH kick ° 


rank ' 


Log{L{0\M)/L,,f) s 


205 


0.1 


2.7 


2.0 


No 


50-50 


0.0 


1 


0.0000000 


229 


0.1 


2.7 


2.0 


Yes 


50-50 


0.0 


2 


-0.057250977 


277 


0.1 


2.7 


2.0 


Yes 


50-50 


0.1 


3 


-0.53945923 


245 


0.1 


2.7 


1.0 


No 


50-50 


0.1 


4 


-1.6570358 


253 


0.1 


2.7 


2.0 


No 


50-50 


0.1 


5 


-1.7356873 


273 


0.1 


2.35 


2.0 


Yes 


50-50 


0.1 


6 


-2.3401947 


269 


0.1 


2.7 


1.0 


Yes 


50-50 


0.1 


12 


-5.6114807 


249 


0.1 


2.35 


2.0 


No 


50-50 


0.1 


10 


-4.1068573 


248 


0.5 


2.7 


1.0 


No 


50-50 


0.1 


55 


-47.292496 


197 


0.1 


2.7 


1.0 


No 


50-50 


0.0 


50 


-42.623398 


241 


0.1 


2.35 


1.0 


No 


50-50 


0.1 


59 


-51.771774 


261 


0.1 


2.7 


0.25 


No 


50-50 


0.1 


81 


-92.800415 


53 


0.1 


2.7 


1.0 


No 


Flat 


0.1 


22 


-15.239967 



^ CE efficiency parameter 

^ Stellar wind strength parameter 

^ 0: CE from Hcrtzsprung gap donor allowed, 1: not allowed 
^ binary mass ratio distribution. 

^ SN kicks given to direct collapse black holes. 0.0 — no SN kick given, 0.1 — small SN kick given 
^ The rank of the model based on the likelihood value 

^ Log of the ratio of the likelihood of the given model to that of the highest likelihood model 




20 40 60 80 200 220 240 260 

Models 1 -96 Models 1 93-288 



280 



Figure 2. Likelihood values for all of the models used in this work. L,.^f is the highest likelihood value among our 192 models. 

We choose a Jeffreys prioiQ on Xij , 

1 



We choose a flat prior on the model parameters, 9, so that 
the posterior is proportional to the likelihood in Eq. dTT ^ 



(8) 



p{e\d) (xp{d\9). 



(11) 



"whence 



2 



Pi^il^) = ^/''^ exp(-A,j) 
Combining Eq. ([9|) and Eq. ([5]), we find that 



pm = n 



'"1,3 



''1,3/ 



1-3 



(1 



■ V. 



!r(i 



(9) 



(10) 



n. 



1,3 ) 



10 



Note that the use of the Jeffreys prior implies that (A; j) = 
+ f • ^ fl3.t prior would have (Aij) = j + 1. Both of 
these priors produce well-defined likelihoods even when j = 
with di^j 7^ 0. The maximum-likelihood estimator, p(Aij|9) = 
S {Xij — Tii^j), while unbiased, produces likelihoods of zero iiriij = 
0, even if only a single count appears in that bin of the data (i.e. 



1). A prior that gives (Xij) 



(i.e. a prior that gives a 



distribution with unbiased mean) is p(Xij) 
in a non-normalizable likelihood when n; , = 



a: 



which results 



4. RESULTS AND DISCUSSION 

Figure 2 plots the likelihood values for each model 
used in this study and shows that only a few models 
are able to accurately reproduce the observed XLFs, so 
this comparison is very useful for eliminating regions of 
our model parameter space. Table 2 lists the six models 
with Log{L{0\M) / Lref) > —3, where Lre/ is the high- 
est likelihood value among our 192 models. These models 
are 205, 229, 277, 245, 253, and 273. These models ah 
have low CE efficiencies {acE = 0.1), a 50-50 mass ratio 
distribution, an IMF exponent of —2.7 (with the excep- 
tion of model 273), and Tywind = 1.0 — 2.0. Recall that our 
CE efficiency parameter really represents acE x A, so low 
values of acE could alternatively be interpreted as these 
systems having a high envelope binding en ergy, which has 
been found to be true for massive stars (jPominik et al.l 
l2012f ). The models that have ?7wind = 2.0, acE = 0.1, 
IMF exponent of —2.7, and a flat q-distribution also have 
fairly high likelihood values (model numbers 13, 37, 61, 



8 



M. Tremmel et al. 



-1 


Jo 0.0< z <0.2 


0.2< z <0.6 


^Siftij. 0.6< z <1.4 


o 
-d. 


.^''in,^ This Paper 


. ^^'"(s. This Paper . 


^^"•^ This Paper J 






F12 


F12 




■4V + T&G08 


"\ + T&G08 


\ + T&G08: 


-3 








-4 
















-5 








-6 


All Galaxies \ 








39 40 41 42 


39 40 41 42 


39 40 41 42 



Log{ Lx (ergs s ) [0.5-2.0 keV] ) 



■X. 


"a 

'o 

Q. 



o 



-1 








-1 


0.0< z <0.4 


^ 0.4<z<1.4 




-2 


This Paper 


This Paper 


X 
03 
T3 


-2 


F12 






•^jt^ + T&G08 


+ T&G08 




-3 






m 

'o 

Q. 


-3 










-4 








-4 














O) 




-5 




o 
—I 


-5 


-6 


Eariy-Type 






-fi 



0.0< z <0.4 


0.4<z<1.4 


.^'^issjj,^ This Paper 


. This Paper 


F12 


^"-s.^ F12 


■ 1 — -N, + T&G08 


>t + T&G08 






Late-Type ^ 





39 40 41 42 

Log( Lx (ergs s"^ 



39 40 41 
[0.5-2.0 keV] ) 



42 



39 40 41 42 

Log( Lx (ergs s"^ 



39 40 41 
[0.5-2.0 keV] ) 



42 



Figure 3. The cyan region shows the area bounded by the six highest UlceUhood models from this study, the grey checkered region shows 
that bounded by the six highest hkehhood models from F12, the black points with error bars are data from T&G08, the sohd black lines 
are from our highest likelihood model (205), and the dashed black lines are from the highest likelihood model from F12 (245). Top: Total 
galaxy population. Bottom left: Early-type galaxies. Bottom right: Late-type galaxies. The XLFs from our highest likelihood models are 
very similar to those from the F12 models, however they underproduce bright early-type galaxies and very bright (Lx > 10*^ ergs s~^) 
late-type galaxies compared with observations. Section 4.1 discusses the causes of these discrepancies in more detail. 



and 85). 

It should be noted that our hkehhood calculation takes 
into account the number of samples in each bin. The 
overall likelihood values are much more sensitive to bins 
with higher sample counts (i.e. those at lower luminos- 

ity). 

Figure 2 and Table 2 show that allowing/not allowing 
CE-HG phases has only a relatively small effect on the 
likelihoods of our best models. In addition, DC BH kicks 
only have an appreciable effect on likelihoods for models 
with lower wind mass loss rates (rewind — 0.25 — 1.0). 
For models with rewind — 2.0, such as those that make 
up the majority of our top models, DC BH kicks have 
little effect. Thus, these two parameters are not very 
well constrained by our analysis. 



4.1. 



Comparison w ith' Fragos et al! [201^) and 
\Tzanavaris & Georgantopoulos (2008) 



F12 use the same PS models used in this work (with the 
inclusion of pure "twins" models) to study the evolution 
of the overall population of XRBs in the Universe. They 
compare with X-ray observations of local galaxies that 
give estimates of the specific X-ray luminosity of XRBs 
in th e local Universe (iLehmer et al]|2010t iBoroson et al.l 
120111: iMineo et al.|[2012[ ). Thev calculate the likelihood 
of each model based on these data and find the six high- 
est likelihood models to be, in order of likelihood, 245, 
229, 269, 205, 249, and 273 (see table 2 for model pa- 
rameters and likelihood values from this study). Figure 
3 compares the six highest likelihood models from this 
work with those of F12. Four out of our top six models 
are also among the six highest likelihood models from 
F12, so it is no surprise that the region bounded by the 



models in this work is very similar to that bounded by 
the models from F12. 

Figure 4 shows XLFs for our highest likelihood model 
(205), with and without hot gas emission, plotted against 
the data from T&G08. Figure 5 plots XLFs similar to 
figure 4, but for the highest likelihood model from F12 
(245). These plots show that our models are able to 
reproduce the redshift evolution of the observed XLFs. 
Consistent with the analysis of T&G08, the XLF evo- 
lution is driven almost entirely by late-type galaxies. 
Our models also reproduce the shape of the early-type 
XLF and the normalization of the late-type XLF, though 
they drastically underproduce bright early-type galax- 
ies and they fail to reproduce the shape of the bright 
(Lx > IQ'^^ergs s'^) end of the late-type XLF. 

Figures 4 and 5 also show that hot gas can have a 
large effect on the shape and normalization of the XLF. 
Further, adding in the hot gas emission suppresses the 
redshift evolution for the early-type galaxy XLF. On 
average, bright early-type galaxies in our models have 
s ~ 60% of their emission come from hot gas while the 
bright late-type galaxies have a 40-50% contribution on 
average. The important role of hot gas emission on the 
XLF means that our simplistic prescriptions for hot gas 
emission add uncertainty to our models and could be a 
major source of the discrepancies between the models 
and observations, particularly for early- type galaxies. 

Another aspect of our models that can account for the 
underproduction of bright early-type galaxies is that our 
PS models only take into account LMXBs formed in the 
field and not those formed dynamically in globular clus- 
ters. Dynamically formed LMXBs are believed to play a 
significant role in old, massive, GC-rich elliptical galax- 
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Figure 4. Luminosity function for the total galaxy population 
using our highest likelihood model, number 205. Solid lines show 
the XLFs from our models with both XRB and hot gas emission. 
The dot-dashed lines show XLFs for just XRB emission. Data 
points and associated error bars are taken from T&G08. The 
model is able to reproduce the evolution and shape of the ob- 
served XLF. Consistent with the analysis of T&G08, the overall 
XLF evolution is driven almost entirely by late- type galaxies. How- 
ever, the model fails to reproduce the correct normalization of the 
Early-type XLF and the shape of the observed XLF of very bright 
(Lx > 10*^ ergs s~^) late-type galaxies. Section 4.1 discusses the 
possible causes for these discrepancies in more detail. Hot gas plays 
an important role in the normalization of the XLF, especially for 
early-type galaxies, where the hot gas contribution also affects the 
XLF evolution with redshift. 

ies. (e.s. iHumphrev fc Buotel[200l iZhang et aLllMl ). 
These LMXB populations can have a significant contri- 
bution to the i ntcKrated X -ray luminosity of bright early- 
type galaxies (|Irwin|[20Q5[ ). Thus, the inclusion of these 
systems could increase the luminosity of our simulated 
early-type galaxies, moving our XLF to the right and 
closer to the observed values. 

For younger, star forming galaxies, our models also 
have only a very basic formula to simulate starburst ac- 
tivity, which can occur, e.g., due to galaxy mergers. This 
would have a significant affect on the HMXB popula- 
tions present in late-type galaxies, and the effect would 
not necessarily be constant with redshift, as merger rates 




Figure 5. Same as figure 4, but for model 245, the highest likeli- 
hood model from F12 and the fifth highest likelihood models from 
this work. These XLFs are similar to that of our highest likelihood 
model, 205. 

may evolve with time (e.g. iConselic^ 120061) . Thus, a 
more detailed SFH is needed to more accurately model 
HMXB populations of late-type galaxies. In addition, 
the higher end of the observed late-type galaxy XLF is 
more at risk from AGN contamination, even with the 
efforts of T&G08 to minimize this effect. A relatively 
small amount of contamination could cause those high 
luminosity data points to move down closer to our model 
predictions. 

4.2. High Redshift Predictions 

Figure 6 plots the X-ray luminosity density from 
normal galaxies as derived from our highest likelihood 
model, 205. The overall evolution (black line in figure 6) 
is very similar to that of the observed SFH of the uni- 
verse. It is also similar to the evolution of the specific 
XRB X-ray luminosity of the universe predicted in F12, 
despite the inclusion here of hot gas emission. This is 
evidence that XRBs drive the overall evolution of the 
normal galaxy X-ray luminosities out to at least z = 4. 
However, our predicted X-ray luminosity density reaches 
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Figure 6. X-ray luminosity density of normal galaxies versus red- 
shift for model 205. The back line shows the contribution from all 
galaxies and the black squares show the total X-ray luminosity den- 
sity from observations presented in T&G08. The X-ray luminosity 
density derived from our models is a factor of 2 — 3 lower than 
the observations. This is due to the fact that our models, relative 
to observations, underproduce bright early-type galaxies and very 
bright (Lx > 10*^) late-type galaxies. Our model follows a simi- 
lar evolution as the SFR density and the X-ray luminosity density 
from XRBs predicted in F12, indicating that XRBs continue to 
drive the evolution of the normal galaxy X-ray emission at higher 
redshifts. However, our models reach a maximum luminosity den- 
sity at 2 ~ 2.5, compared with z ~ 3 for the XRB models in F12. 
This can be attributed to the inclusion of hot gas emission, which 
has been shown to have an appreciable effect the evolution of our 
XLFs when compared with XRB only models. The red, green, and 
blue lines plot the specific X-ray luminosity for different luminosity 
bins. The amount X-ray emission from lower luminosity galaxies 
evolves much less with redshift than that from brighter galaxies. 
In the local Universe, most of the normal galaxy X-ray emission 
comes from lower luminosity galaxies. 

a maximum at z ^ 2.5, which is lower compared with 
the XRB models in F12 that reach a maximum at z ~ 3. 
This can be attributed to the inclusion of hot gas emis- 
sion in our models, which has already been shown to 
have a noticeable effect on the shape and evolution of 
our XLFs, particularly for early-type galaxies. 

Splitting the galaxies into three luminosity bins, 
we find that the evolution of low (10^^ < Lx < 
Wi'^'^ergs s~^) luminosity galaxy emission is small com- 
pared with the evolution for higher luminosity galaxies, 
which varies by an order of magnitude on the range 
2: = to z = 4. The most luminous galaxies (Lx > 
Vii^^ergs s~^) reach a maximum around z = 3, which 
also approximately corresponds to the time of maxi- 
mum SFR density in the Universe. Galaxies in the 
range 10*° < Lx < Vd'^^ergs reach a maximum 
around z — 2. In the local Universe, the low lumi- 
nosity galaxies dominate the normal galaxy X-ray emis- 
sion. We do not go to higher redshift here because our 
hot gas emission prescription relies on galaxy morphol- 
ogy, which becomes h arder to classify at higher redshifts 
(|van den Bergh|[2Q0l . 

4.3. Effects of parameters on XLFs 

Figure 7 shows XLFs for different models compared 
to model 245. Each model is chosen to encapsulate the 
effect that each parameter has on the shape of the XLF. 




39 40 41 

Log( (ergs s"') [0.5-2.0 keV] ) 



Figure 7. the highest likelihood model from F12(245, black line) 
compared with other models to illustrate the effects of different 
parameters on the shape of the XLF. Different parameters have 
varying effects on the shape of the XLFs. Shown here are mod- 
els with higher acE (248), CE-HG allowed (269), no BH kicks 
(197), steeper IMF (241), lower wind mass loss (261), and a fiat 
distribution (53). 

Several parameters have significant effects on the shape 
of the XLFs, while others have only minimal effects. 

The common envelope efficiency parameter (acs) dic- 
tates how efficiently orbital energy is converted to ther- 
mal energy that will expel the envelope. A lower effi- 
ciency means that it will take more orbital energy to ex- 
pel the envelope. This parameter mainly effects LMXBs, 
as most LMXBs formed in the field must go through a CE 
phase. The CE phase plays an important role in making 
the orbit close enough to allow for RLO, but a lower CE 
efficiency leads to even more orbital decay and a higher 
rate of mergers, overall decreasing the rate of LMXB 
formation. This effect can be seen by comparing models 
245 and 248. HMXBs are not as strongly affected by 
changes in acEi as they have other forma tion channels 
available that do not invo lve CE phases ([Linden et al.l 
l201Clt iValsecchi et al.l 120101 ) . This parameter mainly af- 
fects the lower luminosity end of the XLF, where a higher 
CE efficiency increases the number of bright galaxies due 
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to an increased LMXB population. 

Wind mass loss rates effect the evolution of high mass 
stars in two major competing ways. Higher wind mass 
loss rates will increase the accretion rates of wind-fed 
HMXBs, increasing their luminosity. On the contrary, 
lower winds will result in a lower overall mass loss of 
the primary star and hence increase the formation rate 
of massive BHs. BH-XRBs tend to be more luminous 
than XRBs with NS accretors. This is because, on the 
one hand, they can form stable RLO XRB systems with 
massive companions and, on the other hand, BHs drive 
higher accretion rates due to their higher mass and, 
therefore, higher Eddington limits. In this way, weaker 
stellar winds can increase the luminosities of both LMXB 
and HMXB populations. Comparing models 245 and 
261, we see that weaker stellar winds increases the num- 
ber of bright early and late-type galaxies, so the latter 
effect is dominant. 

Changing the initial binary mass ratio between a flat 
distribution and a 50% twins and 50% fiat distribution ef- 
fects both the HMXB and LMXB populations. As stated 
earlier, the binary mass ratio effects the secondary star 
in the system, which will eventually become the donor 
star in most cases. All XRBs are accreting mass onto a 
compact object, which can only come from a high mass 
progenitor. Most LMXBs require a high initial mass ra- 
tio in order to ensure a high mass primary star that will 
evolve into a BH or NS with a lower mass companion. 
The "50-50" distribution adopted here forces mass ratios 
close to 1 for half the binary population. This will de- 
crease the LMXB population while increasing the HMXB 
population, as HMXBs require ratios closer to 1. Chang- 
ing from "50-50" to a flat distribution (comparing mod- 
els 245 and 53) increases the lower luminosity end of the 
early and late-type galaxy XLFs and slightly decreases 
the number of very bright late-type galaxies galaxies. 

Allowing small natal kicks for direct collapse BHs af- 
fects the HMXB and LMXB population in two compet- 
ing ways. On the one hand, nat al kicks can enh ance the 
formation of RLO-HMXBs (Li nden et al.ll2010t ). On the 
other hand, natal kicks inject energy into the binary sys- 
tem and could result in the widening or the complete 
disruption of the system, thereby decreasi ng the forma- 
tion H MXB and LMXB with BH accretors. iLinden et all 
(|2009[ ) find that imparting small natal kicks to DC BHs 
is necessary in order to reproduce the lack of observed 
wide orbit BH-XRBs. Comparing models 245 and 197, 
not allowing natal kicks increases the high luminosity end 
of the late-type galaxy XLF and has little effect on the 
early-type galaxy population. 

Within our grid of PS models, we also have the 
IMF power law exponent as a free parameter, allow- 
ing it to be either —2.35 or —2.7. It is instructive to 
note that the IMF referred t o in this work represents 
the integrated galaxy I MF (jWeidner fc Kroupal 120051: 
iKroupa fc Weidneill2003[ ) and that this work only probes 
the high mass region of the IMF, since the primary stars 
that are created via sampling the IMF must be massive 
enough to form a BH or NS. The slope of the power law 
at the high-mass end of the IMF affects the population 
of XRBs in a way similar to stellar winds, in the sense 
that a flatter IMF will produce relatively more massive 
BHs compared to a steeper one. A flatter IMF will result 
in more bright LMXBs and HMXBs. This affect can be 



seen by comparing models 245 and 241. As expected, a 
flatter IMF results in a higher number of bright galaxies, 
though the very high luminosity end of the early-type 
galaxy XLF is not strongly affected by this parameter. 
This is not surprising, as figure 4 shows that hot gas 
emission dominates the high luminosity end of the XLF 
for early- type galaxies. 

Finally, we found that allowing or not for all pos- 
sible outcomes in CE phases with donor stars in the 
Hertzsprung gap has only a small affect on the galaxy 
XLFs, slightly increasing the number of bright late- type 
galaxies (compare models 245 and 269). Thus, although 
this parameter affec ts the shape of t he XLF of individual 
XRBs in a galaxy ()Luo et al.ll2012[ ). it has a negligible 
effect on the integrated X-ray luminosity of a galaxy. 

The models that we find to agree best with observa- 
tions have direct collapse BH natal kicks, low CE ef- 
ficiency, steep IMFs, -qwind — 1-0 — 2.0, and a 50-50 
mass ratio distribution. As outlined above, these pa- 
rameters make for a limited population of both LMXBs 
and BH-XRBs. The HMXB population is limited by the 
steeper IMF, BH natal kicks, and higher winds, but also 
strengthened by the 50-50 mass ratio distribution. How- 
ever, the latter effect on the XLF is much weaker, as 
shown in figure 7. 

5. CONCLUSIONS 

Using data from the Millennium C'osmological Simula- 
tion and the semi-analytical analysis conducted by Gil 
in tandem with the binary population synthesis code, 
StarTrack, we simulated the population of XRBs within 
normal galaxies in a large volume of the Universe from 
z = to ~ 20. Assuming that galaxy X-ray emission is 
solely due to XRBs and hot gas, we calculated the in- 
tegrated X-ray luminosity of each galaxy in this cosmic 
volume and compared the resulting galaxy XLFs to the 
observational XLFs of T&G08. 

In this paper, we presented data from 192 binary pop- 
ulation synthesis models, varying parameters that have 
the largest effect on binary star evolution (see Table 1 
for a list of the parameters) . We use a likelihood calcu- 
lation method to compare each model with the results 
from T&G08. From this analysis we find that only a 
few of our models are able to accurately reproduce the 
most recent observations of X-ray bright normal galaxies. 
Our highest likelihood models are also among the highest 
likelihood models from a separate analysis presented in 
F12. This confirms that our results are consistent with 
their separate analysis, which compares these same mod- 
els with the observed overall emission from LMXB and 
HMXB populations in the local Universe. To have only 
^ 10 models from our 192 model grid that accurately 
match observations in two separate analysis shows that 
we are able to provide remarkable and self consistent 
constrains on the XRB parameter space. 

We find that our highest likelihood models are those 
with a lower LMXB population due to a low CE efficiency 
and a 50-50 mass ratio distribution, and a lower BH- 
XRB and HMXB population due to higher winds, and 
a steeper IMF. Our models do well in reproducing the 
normalization and evolution of the total and late-type 
galaxy XLFs, as well as the evolution and shape of the 
early-type XLF. 

Our models show that hot gas emission has a large 
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effect on the shape of the XLFs, and it signifieantly af- 
fects the redshift evolution of the early-type galaxy XLF, 
causing it to remain nearly constant out to z = 1.4. 

We show that the observed redshift evolution of the 
normal galaxy XLF continues out to higher redshift, with 
the specific normal galaxy X-ray luminosity evolving in 
a way similar to the SFH of the imivcrse and consistent 
with the evolution of XRB emission found in F12. This 
is evidence that the XLF evolution is driven by XRB 
evolution even out to higher redshifts. Our models also 
show that hot gas emission causes the point of maximum 
normal galaxy X-ray luminosity density to shift to lower 
redshift compared with the XRB models in F12. 

However, dcispitc these many successes, our models do 
not perfectly reproduce the observed XLFs. In partic- 
ular, they fail to reproduce the observed normalization 
of the early-type galaxy XLF, greatly underestimating 
the number of bright early-type galaxies. Our highest 
likelihood models also fail to reproduce the shape of the 
high (Lx > lO'^^ergs s~^) luminosity end of the late-type 
galaxy XLF, particularly for higher redshifts. 

Our models have limitations that may have caused 
these discrepancies. For one, we do not take into ac- 
count dynamically formed LMXBs, which could signifi- 
cantly increase the normalization of the model early-type 
galaxy XLF. For late-type galaxies, the XRB luminosi- 
ties have a higher contribution from HMXB populations, 
which are very sensitive to evolving SFHs. However, the 
SFHs used from the Gil catalog are limited in their de- 
tail and our method for simulating the effect of starbursts 
is very rudimentary. In addition, our prescription for 
hot gas, though based on observations, is very basic and 
could add inaccuracy to our X-ray luminosities. A more 
detailed model is needed to more accurately model the 
hot gas emission, particularly in early-type galaxies. 

In addition to limitations in our models, the observa- 
tions of very bright galaxies are subject to the possibility 
of AGN contamination, which could artificially increase 
the observed high luminosity data points from T&G08. 

Despite these shortcomings, this work represents a 
first careful attempt to study how XRBs control the Lx 
distributions of different types of galaxies. As such it, 
provides an important theoretical base for future X-ray 
observations of normal galaxies at high redshift. It also 
shows that XRB populations are closely linked with 
the growth of galaxies. This work lays the ground for 
future work using X-ray observations and cosmological 
simulations of galaxies to provide a new way to constrain 
our models of binary evolution, as well as study the 
role played by XRBs in galaxy formation and evolution 
through feedback processes. 
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Table 3 

A complete list of all PS models used in this work. 



Model 




IMF exponent 


rj ^ 
Wwind 




q distribution " 


T\*^ T>TT 1 ■ 1 P 

UC BH kick 


1 f 
rank 


Log(L{0\M)IL^^^) s 


1 


0.1 


2.35 


1.0 


No 


Flat 


0.0 


105 


-172.65289 


2 


0.2 


2.35 


1.0 


No 


Flat 


0.0 


127 


-298.51453 


3 


0.3 


2.35 


1.0 


No 


Flat 


0.0 


142 


-439.21210 


4 


0.5 


2.35 


1.0 


No 


Flat 


0.0 


161 


-727.54431 


5 


0.1 


2.7 


1.0 


No 


Flat 


0.0 


42 


-33.895195 


6 


0.2 


2.7 


1.0 


No 


Flat 


0.0 


72 


-75.746765 


7 


0.3 


2.7 


1.0 


No 


Flat 


0.0 


89 


-111.71681 


8 


0.5 


2.7 


1.0 


No 


Flat 


0.0 


106 


-176.81656 


9 


0.1 


2.35 


2.0 


No 


Flat 


0.0 


29 


-16.885155 


10 


0.2 


2.35 


2.0 


No 


Flat 


0.0 


69 


-72.331543 


11 


0.3 


2.35 


2.0 


No 


Flat 


0.0 


93 


-129.83803 


12 


0.5 


2.35 


2.0 


No 


Flat 


0.0 


90 


-113.39149 


13 


0.1 


2.7 


2.0 


No 


Flat 


0.0 


11 


-4.6178589 


14 


0.2 


2.7 


2.0 


No 


Flat 


0.0 


37 


-26.552711 


15 


0.3 


2.7 


2.0 


No 


Flat 


0.0 


49 


-41.962570 


16 


0.5 


2.7 


2.0 


No 


Flat 


0.0 


44 


-36.056236 


17 


0.1 


2.35 


0.25 


No 


Flat 


0.0 


165 


-779.69275 


18 


0.2 


2.35 


0.25 


No 


Flat 


0.0 


176 


-1189.3790 


19 


0.3 


2.35 


0.25 


No 


Flat 


0.0 


184 


-1683.3291 


20 


0.5 


2.35 


0.25 


No 


Flat 


0.0 


190 


-2598.7773 


21 


0.1 


2.7 


0.25 


No 


Flat 


0.0 


119 


-247.78528 


22 


0.2 


2.7 


0.25 


No 


Flat 


0.0 


141 


-423.04449 


23 


0.3 


2.7 


0.25 


No 


Flat 


0.0 


152 


-572.43054 


24 


0.5 


2.7 


0.25 


No 


Flat 


0.0 


169 


-904.26648 


25 


0.1 


2.35 


1.0 


Yes 


Flat 


0.0 


112 


-204.92599 


26 


0.2 


2.35 


1.0 


Yes 


Flat 


0.0 


136 


-375.42163 


27 


0.3 


2.35 


1.0 


Yes 


Flat 


0.0 


151 


-568.51801 


28 


0.5 


2.35 


1.0 


Yes 


Flat 


0.0 


173 


-1074.5802 


29 


0.1 


2.7 


1.0 


Yes 


Flat 


0.0 


48 


-40.053345 


30 


0.2 


2.7 


1.0 


Yes 


Flat 


0.0 


82 


-94.520035 


31 


0.3 


2.7 


1.0 


Yes 


Flat 


0.0 


101 


-160.10526 


32 


0.5 


2.7 


1.0 


Yes 


Flat 


0.0 


125 


-288.00360 


33 


0.1 


2.35 


2.0 


Yes 


Flat 


0.0 


32 


-19.284058 


34 


0.2 


2.35 


2.0 


Yes 


Flat 


0.0 


80 


-92.763794 


35 


0.3 


2.35 


2.0 


Yes 


Flat 


0.0 


100 


-156.03748 


36 


0.5 


2.35 


2.0 


Yes 


Flat 


0.0 


104 


-169.66629 


37 


0.1 


2.7 


2.0 


Yes 


Flat 


0.0 


8 


-3.8255920 


38 


0.2 


2.7 


2.0 


Yes 


Flat 


0.0 


38 


-26.668427 


39 


0.3 


2.7 


2.0 


Yes 


Flat 


0.0 


57 


-48.150604 


40 


0.5 


2.7 


2.0 


Yes 


Flat 


0.0 


54 


-46.984512 


41 


0.1 


2.35 


0.25 


Yes 


Flat 


0.0 


164 


-776.80261 


42 


0.2 


2.35 


0.25 


Yes 


Flat 


0.0 


179 


-1404.6885 


43 


0.3 


2.35 


0.25 


Yes 


Flat 


0.0 


186 


-2075.5227 


44 


0.5 


2.35 


0.25 


Yes 


Flat 


0.0 


192 


-3790.4802 


45 


0.1 


2.7 


0.25 


Yes 


Flat 


0.0 


122 


-270.38574 


46 


0.2 


2.7 


0.25 


Yes 


Flat 


0.0 


147 


-506.66522 


47 


0.3 


2.7 


0.25 


Yes 


Flat 


0.0 


162 


-760.54401 


48 


0.5 


2.7 


0.25 


Yes 


Flat 


0.0 


182 


-1476.1992 


49 


0.1 


2.35 


1.0 


No 


Flat 


0.1 


74 


-78.037476 


50 


0.2 


2.35 


1.0 


No 


Flat 


0.1 


108 


-191.24628 


51 


0.3 


2.35 


1.0 


No 


Flat 


0.1 


132 


-337.98636 


52 


0.5 


2.35 


1.0 


No 


Flat 


0.1 


154 


-637.43152 


53 


0.1 


2.7 


1.0 


No 


Flat 


0.1 


22 


-15.239967 


54 


0.2 


2.7 


1.0 


No 


Flat 


0.1 


58 


-48.793671 


55 


0.3 


2.7 


1.0 


No 


Flat 


0.1 


78 


-87.307587 


56 


0.5 


2.7 


1.0 


No 


Flat 


0.1 


94 


-133.84480 


57 


0.1 


2.35 


2.0 


No 


Flat 


0.1 


30 


-17.080589 


58 


0.2 


2.35 


2.0 


No 


Flat 


0.1 


68 


-69.020309 


59 


0.3 


2.35 


2.0 


No 


Flat 


0.1 


91 


-117.72034 


60 


0.5 


2.35 


2.0 


No 


Flat 


0.1 


86 


-104.91006 


61 


0.1 


2.7 


2.0 


No 


Flat 


0.1 


7 


-3.7329865 


62 


0.2 


2.7 


2.0 


No 


Flat 


0.1 


34 


-22.856674 


63 


0.3 


2.7 


2.0 


No 


Flat 


0.1 


46 


-39.461548 


64 


0.5 


2.7 


2.0 


No 


Flat 


0.1 


43 


-35.445259 


65 


0.1 


2.35 


0.25 


No 


Flat 


0.1 


155 


-646.17456 


66 


0.2 


2.35 


0.25 


No 


Flat 


0.1 


174 


-1132.3506 


67 


0.3 


2.35 


0.25 


No 


Flat 


0.1 


183 


-1587.7677 


68 


0.5 


2.35 


0.25 


No 


Flat 


0.1 


189 


-2579.9290 


69 


0.1 


2.7 


0.25 


No 


Flat 


0.1 


103 


-164.22797 


70 


0.2 


2.7 


0.25 


No 


Flat 


0.1 


131 


-330.11389 


71 


0.3 


2.7 


0.25 


No 


Flat 


0.1 


145 


-480.67297 


72 


0.5 


2.7 


0.25 


No 


Flat 


0.1 


168 


-889.29907 


73 


0.1 


2.35 


1.0 


Yes 


Flat 


0.1 


88 


-108.42148 


74 


0.2 


2.35 


1.0 


Yes 


Flat 


0.1 


123 


-275.61371 


75 


0.3 


2.35 


1.0 


Yes 


Flat 


0.1 


146 


-500.68408 
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a 


AiVlJT tJAULIIltJIllj 


n ■ 1 b 
'/wind 


CE-HG " 


vj^ UloLllUUtlUli 


DC, RH VirV " 


ran 


T,nn( T,(n\M\ 1 T, 


76 


0.5 


2.35 


1.0 


Yes 


Flat 


0.1 


171 


-1027.9811 


77 


0.1 


2.7 


1.0 


Yes 


Flat 


0.1 


27 


-16.603905 


78 


0.2 


2.7 


1.0 


Yes 


Flat 


0.1 


67 


-68.651489 


79 


0.3 


2.7 


1.0 


Yes 


Flat 


0.1 


92 


-126.54276 


80 


0.5 


2.7 


1.0 


Yes 


Flat 


0.1 


117 


-236.25699 


81 


0.1 


2.35 


2.0 


Yes 


Flat 


0.1 


20 


-13.885048 


82 


0.2 


2.35 


2.0 


Yes 


Flat 


0.1 


75 


-80.561569 


83 


0.3 


2.35 


2.0 


Yes 


Flat 


0.1 


95 


-136.36365 


84 


0.5 


2.35 


2.0 


Yes 


Flat 


0.1 


102 


-161.30212 


85 


0.1 


2.7 


2.0 


Yes 


Flat 


0.1 


9 


-3.8542175 


86 


0.2 


2.7 


2.0 


Yes 


Flat 


0.1 


35 


-25.809174 


87 


0.3 


2.7 


2.0 


Yes 


Flat 


0.1 


53 


-45.016434 


88 


0.5 


2.7 


2.0 


Yes 


Flat 


0.1 


52 


-44.284752 


89 


0.1 


2.35 


0.25 


Yes 


Flat 


0.1 


156 


-649.58960 


90 


0.2 


2.35 


0.25 


Yes 


Flat 


0.1 


177 


-1272.5421 


91 


0.3 


2.35 


0.25 


Yes 


Flat 


0.1 


185 


-2018.7471 


92 


0.5 


2.35 


0.25 


Yes 


Flat 


0.1 


191 


-3722.3474 


93 


0.1 


2.7 


0.25 


Yes 


Flat 


0.1 


107 


-184.78189 


94 


0.2 


2.7 


0.25 


Yes 


Flat 


0.1 


139 


-405.31018 


95 


0.3 


2.7 


0.25 


Yes 


Flat 


0.1 


159 


-687.57117 


96 


0.5 


2.7 


0.25 


Yes 


Flat 


0.1 


180 


-1430.8417 


193 


0.1 


2.35 


1.0 


No 


50-50 


0.0 


109 


-191.83087 


194 


0.2 


2.35 


1.0 


No 


50-50 


0.0 


118 


-243.62000 


195 


0.3 


2.35 


1.0 


No 


50-50 


0.0 


126 


-288.46088 


196 


0.5 


2.35 


1.0 


No 


50-50 


0.0 


138 


-397.25439 


197 


0.1 


2.7 


1.0 


No 


50-50 


0.0 


50 


-42.623398 


198 


0.2 


2.7 


1.0 


No 


50-50 


0.0 


65 


-60.471497 


199 


0.3 


2.7 


1.0 


No 


50-50 


0.0 


70 


-73.344910 


200 


0.5 


2.7 


1.0 


No 


50-50 


0.0 


84 


-97.657059 


201 


0.1 


2.35 


2.0 


No 


50-50 


0.0 


19 


-13.511047 


202 


0.2 


2.35 


2.0 


No 


50-50 


0.0 


45 


-37.678787 


203 


0.3 


2.35 


2.0 


No 


50-50 


0.0 


64 


-60.211609 


204 


0.5 


2.35 


2.0 


No 


50-50 


0.0 


63 


-54.762085 


205 


0.1 


2.7 


2.0 


No 


50-50 


0.0 


1 


0.0000000 


206 


0.2 


2.7 


2.0 


No 


50-50 


0.0 


14 


-8.9467392 


207 


0.3 


2.7 


2.0 


No 


50-50 


0.0 


25 


-16.300880 


208 


0.5 


2.7 


2.0 


No 


50-50 


0.0 


18 


-12.898163 


209 


0.1 


2.35 


0.25 


No 


50-50 


0.0 


149 


-523.74426 


210 


0.2 


2.35 


0.25 


No 


50-50 


0.0 


160 


-700.39764 


211 


0.3 


2.35 


0.25 


No 


50-50 


0.0 


170 


-943.16479 


212 


0.5 


2.35 


0.25 


No 


50-50 


0.0 


181 


-1438.5382 


213 


0.1 


2.7 


0.25 


No 


50-50 


0.0 


111 


-198.82269 


214 


0.2 


2.7 


0.25 


No 


50-50 


0.0 


121 


-264.05933 


215 


0.3 


2.7 


0.25 


No 


50-50 


0.0 


130 


-328.01355 


216 


0.5 


2.7 


0.25 


No 


50-50 


0.0 


144 


-471.44672 


217 


0.1 


2.35 


1.0 


Yes 


50-50 


0.0 


115 


-217.41257 


218 


0.2 


2.35 


1.0 


Yes 


50-50 


0.0 


124 


-279.08957 


219 


0.3 


2.35 


1.0 


Yes 


50-50 


0.0 


133 


-351.58136 


220 


0.5 


2.35 


1.0 


Yes 


50-50 


0.0 


150 


-558.37683 


221 


0.1 


2.7 


1.0 


Yes 


50-50 


0.0 


60 


-52.107895 


222 


0.2 


2.7 


1.0 


Yes 


50-50 


0.0 


71 


-74.032425 


223 


0.3 


2.7 


1.0 


Yes 


50-50 


0.0 


85 


-99.106094 


224 


0.5 


2.7 


1.0 


Yes 


50-50 


0.0 


96 


-141.66830 


225 


0.1 


2.35 


2.0 


Yes 


50-50 


0.0 


33 


-22.380219 


226 


0.2 


2.35 


2.0 


Yes 


50-50 


0.0 


62 


-53.372421 


227 


0.3 


2.35 


2.0 


Yes 


50-50 


0.0 


73 


-77.036163 


228 


0.5 


2.35 


2.0 


Yes 


50-50 


0.0 


77 


-85.713593 


229 


0.1 


2.7 


2.0 


Yes 


50-50 


0.0 


2 


-0.057250977 


230 


0.2 


2.7 


2.0 


Yes 


50-50 


0.0 


16 


-9.4789734 


231 


0.3 


2.7 


2.0 


Yes 


50-50 


0.0 


31 


-17.890610 


232 


0.5 


2.7 


2.0 


Yes 


50-50 


0.0 


28 


-16.608932 


233 


0.1 


2.35 


0.25 


Yes 


50-50 


0.0 


148 


-510.56641 


234 


0.2 


2.35 


0.25 


Yes 


50-50 


0.0 


167 


-813.52740 


235 


0.3 


2.35 


0.25 


Yes 


50-50 


0.0 


175 


-1156.5618 


236 


0.5 


2.35 


0.25 


Yes 


50-50 


0.0 


188 


-2207.2117 


237 


0.1 


2.7 


0.25 


Yes 


50-50 


0.0 


113 


-206.70996 


238 


0.2 


2.7 


0.25 


Yes 


50-50 


0.0 


128 


-299.82571 


239 


0.3 


2.7 


0.25 


Yes 


50-50 


0.0 


140 


-417.71179 


240 


0.5 


2.7 


0.25 


Yes 


50-50 


0.0 


163 


-771.67303 


241 


0.1 


2.35 


1.0 


No 


50-50 


0.1 


59 


-51.771774 


242 


0.2 


2.35 


1.0 


No 


50-50 


0.1 


83 


-95.560333 


243 


0.3 


2.35 


1.0 


No 


50-50 


0.1 


97 


-147.48036 


244 


0.5 


2.35 


1.0 


No 


50-50 


0.1 


120 


-257.68997 


245 


0.1 


2.7 


1.0 


No 


50-50 


0.1 


4 


-1.6570358 


246 


0.2 


2.7 


1.0 


No 


50-50 


0.1 


23 


-15.514267 


247 


0.3 


2.7 


1.0 


No 


50-50 


0.1 


40 


-29.969849 
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Model acE " IMF exponent rj^ind ^ CE-HG " q distribution <1 DC BH kick " rank ^ Log(L(0|M)/L^e/) ^ 



248 


0.5 


2.7 


1.0 


No 


50-50 


0.1 


55 


47 9Q9ZLC1R 


249 


0.1 


2.35 


2.0 


No 


50-50 


0.1 


10 


I.IUUOO t o 


250 


0.2 


2.35 


2.0 


No 


50-50 


0.1 


36 




251 


0.3 


2.35 


2.0 


No 


50-50 


0.1 


51 


-^O.UOOO ( '-t 


252 


0.5 


2.35 


2.0 


No 


50-50 


0.1 


47 


oy .0U100U 


253 


0.1 


2.7 


2.0 


No 


50-50 


0.1 


5 


- 1 . ( OOUO ( o 


254 


0.2 


2.7 


2.0 


No 


50-50 


0.1 


13 


f .uououuo 


ZOO 


u.o 


9 7 


9 n 
z.u 


No 


ou-ou 


n 1 

U. 1 


21 


-l^.OZOZOO 


ZOO 


u.o 


9 7 


9 n 
z.u 


IMO 


ou-ou 


n 1 

U. 1 


1 7 
1 ( 


-ii.yyZ'rioo 


257 


0.1 


2.35 


0.25 


No 


50-50 


0.1 


135 


-•^fifi Q9471 


ZOO 


n 9 


9 'i^ 


n 9^^ 

u. zo 


No 


ou-ou 


n 1 

U. 1 


lOO 






u.o 


z.oo 


u. zo 


No 


ou-ou 


U. 1 


lUU 


-oU'y^.oui^o 




u.o 


z.oo 


u. zo 


No 


ou-ou 


U. 1 


1 7S 
1 i o 


-lO'y: 1 .yzoi 


261 


0.1 


2.7 


0.25 


No 


50-50 


0.1 


81 


Q9 ^0041 


262 


n 9 
u.z 


9 7 
Z. i 


u.zo 


No 


ou-ou 


n 1 

U. 1 


QR 
yo 


1 f^*^ 9Q4'^'^ 
-loo.zy'ioo 


9fi'^ 

ZlJO 


u.o 


9 7 

Z. J 


n 9^^ 

u.zo 


No 


ou-ou 


1 


114 


-ziz.oooyu 


264 


0.5 


2.7 


0.25 


No 


50-50 


0.1 


137 


-oyz.uooyo 


265 


0.1 


2.35 


1.0 


Yes 


50-50 


0.1 


79 


-oo.youuou 


9fifi 
ZOU 


n 9 
u.z 


z.oo 


1 n 

l.U 


1 cb 


ou-ou 


n 1 

U. 1 


QQ 


-lOO. / U / O ( 


267 


0.3 


2.35 


1.0 


Yes 


50-50 


0.1 


116 


-ZOO.Z 1 uou 


268 


0.5 


2.35 


1.0 


Yes 


50-50 


0.1 


143 


-447 R7(^f,^ 

I .O 1 UUl 




n 1 

U. J. 


9 7 

Z. 1 


1 n 

l.U 


I CO 


ou-ou 


n 1 

U. 1 


12 


1 4Qn7 

-O.Oll^OU / 


97n 


n 9 

u.z 


9 7 
z. r 


1 
l.U 


I CO 


ou-ou 


1 
U. 1 


oy 


97 4'^9fiS9 

-Z / .^OZOoZ 


271 


0.3 


2.7 


1.0 


Yes 


50-50 


0.1 


56 


47 QS94 


272 


u.o 


9 7 

Z. f 


1 n 

l.U 


Yes 


ou-ou 


U. 1 


7fi 


-oo. you'ioz 




n 1 

U. J- 


z.oo 


9 n 

z.u 


Yes 


ou-ou 


n 1 

U. 1 


u 


9 "^401 Q47 

-Z.04:Uiy4: I 


274 


0.2 


2.35 


2.0 


Yes 


50-50 


0.1 


41 


'il 'iQ9n'^fi 
-o 1 .oy zuou 


275 


0.3 


2.35 


2.0 


Yes 


50-50 


0.1 


61 


074701 
-oz . y 1 4: 1 u 1 


97R 


u.o 


z.oo 


9 n 
z.u 


Yes 


ou-ou 


U. 1 


uu 


-uo.oyuulo 


977 


n 1 

U. -L 


9 7 

Z. i 


9 n 
z.u 


Yes 


ou-ou 


U. 1 


o 




278 


0.2 


2.7 


2.0 


Yes 


50-50 


0.1 


15 


-9.2003784 


279 


0.3 


2.7 


2.0 


Yes 


50-50 


0.1 


24 


-16.015152 


280 


0.5 


2.7 


2.0 


Yes 


50-50 


0.1 


26 


-16.512665 


281 


0.1 


2.35 


0.25 


Yes 


50-50 


0.1 


134 


-364.58423 


282 


0.2 


2.35 


0.25 


Yes 


50-50 


0.1 


157 


-666.35565 


283 


0.3 


2.35 


0.25 


Yes 


50-50 


0.1 


172 


-1048.5062 


284 


0.5 


2.35 


0.25 


Yes 


50-50 


0.1 


187 


-2083.5793 


285 


0.1 


2.7 


0.25 


Yes 


50-50 


0.1 


87 


-105.98897 


286 


0.2 


2.7 


0.25 


Yes 


50-50 


0.1 


110 


-193.32385 


287 


0.3 


2.7 


0.25 


Yes 


50-50 


0.1 


129 


-315.21362 


288 


0.5 


2.7 


0.25 


Yes 


50-50 


0.1 


158 


-687.04688 



^ CE efficiency parameter 

^ Stellar wind strength parameter 

^ 0: CE from Hcrtzsprung gap donor allowed, 1: not allowed 
^ binary mass ratio distribution. 

^ SN kicks given to direct collapse black holes. 0.0 — no SN kick given, 0.1 = small SN kick given 
f 

The rank of the model based on the likelihood value 
^ Log of the ratio of the likelihood of the given model to that of the highest likelihood model 



